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Massless Dirac particles on the helicoid are theoretically investigated. With its possible application being 
helical graphene, we explore how the peculiarities of Dirac particles appear on the curved, screw-symmetric 
surface. Zweibein is used to derive the massless Dirac equation on the helicoid and on general curved surfaces. 
We show that bound states of massless Dirac electrons on the helicoid are shown to be absent, and thus the 
system is fully characterized by the scattering probabilities and the phase shifts. We obtained these quantities 
from numerically calculated wavefunctions. We find the local density of states and the phase shifts behave 
characteristically around the axis of the helicoid. Bound states of massive Dirac electrons on the surface are 
also shown to be absent as an extension of the above result on massless Dirac electrons. A comparison with the 
non-relativistic case is also made. 

PACS numbers: 73.20.At, 02.40.-k 


I. INTRODUCTION 


Graphene and massless Dirac fermions on it have gained 
persistent attention over the last few decades as a way of look¬ 
ing at relativity through condensed matter physics 1 ' While 
there are a number of studies on graphene, only a few focus on 
curved graphene and its electronic stuructures. It is, therefore, 
worth asking the following question: how does the electronic 
structure of graphene change if we deform it into a curved 
surface? This is interesting in a number of ways, (i) Non- 
relativistic electrons on periodic curved surfaces have distinc¬ 
tive electronic structures? and we expect something even more 
interesting with Dirac fermions due to their spinor wavefunc- 
tions rather than scalar ones, (ii) It is also intriguing to see 
how the presence of positive and negative energy states of 
Dirac fermions exerts its effect on curved surfaces. In other 
words how the Klein paradox 3 } i.e., potential barriers cannot 
reflect massless Dirac particles, takes its form on curved sur¬ 
faces. 


These have motivated us to look into Dirac particles on the 
helicoid, one of the simplest curved surfaces periodic in one 
direction's, xhe surface is also simple in that it is a minimal 
surface, i.e., its mean curvature vanishes everywhere - it is the 
only ruled minimal surface other than the plane. Possible ap¬ 
plications to physical systems include graphite with a screw 
dislocation. Horn for example reported a spiral growth pat¬ 
terns on natural graphite as early as 1952 5 , followed by more 
recent studies by Rakovan and JaszczakP. There is also a the¬ 
oretical study by Bird and Preston of a spiral graphite in terms 
of Berry’s phastP. 


The rest of the paper is organized as follows. In Section [Ii] 
we formulate the Dirac equation on general curved surfaces, 
and then in Section III apply the formalism to the helicoid. 
In Section IIIICI we show that there are no bound states of 


massless Dirac fermions on the helicoid and discuss the scat¬ 
tering amplitudes and the phase shifts of an electron off the 
spiral axis in terms of partial waves. We go on to discuss 
the local density of states in Section[IV] In Section[V]we ex¬ 
amine the (non-)existence of bound states of massive Dirac 
fermions on the helicoid. In Appendix[B]comparison with the 
non-relativistic case is briefly made. 


H. DIRAC EQUATION ON CURVED SURFACES 

The massless Dirac equation on a flat two-dimensional sur¬ 
face. 


= Ei/s, (1) 

describes the low-energy behaviour of electrons on graphene 1 . 
Here cr^ (p = 1,2) are the Pauli matrices, t ft is the spinor 
wavefunction, and E is the eigenenergjC 12 . Here and hereafter 
repeated indices are implicitely summed over. We, following 
Birrell and Davies} generalize the Dirac equation to curved 
surfaces by replacing derivatives with covariant derivatives, 
and by introducing zweibeins, e^: 

icr 0 e b D = Etp (a , p = 1, 2). (2) 

We take a convention in which the metric is given by g /IV = 
efevOab with rjab = diag(+, -, -) being the Minkowskian 
metric. The covariant derivatives for spinors are given by 
D fl - if, + T /( , where 

r M = \s ab e v a g pv D l ,e P b and S ab = ^[cr a , cr b ]. (3) 
The action of covariant derivatives on zweibeins is 


D v e; = d v e; + e b ojl 




A 

v/r 


(4) 
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Here u> a b are the spin connection coefficients defined via 

d(e“dx?) = -<jf b A ( e b v dx v ), (5) 


where the coordinate v spirals about the axis of the helicoid, 
while u shoots out normal to v from the axis with a being the 
pitch of the spiral. The metric in this parametrization is 


while are the ChristofFel symbols. We require the spin 
connection coefficients to be antisymmetric in a and b so as to 
avoid the arbitrariness in the above definition. 

Every surface has global isothermal coordinates, i.e., every 
metric on any surfaces can be represented in a diagonal form 
upon a suitable coordinate transformation 9 ^ 


The zweibeins and ChristofFel symbols then read, in matrix 
forms. 


e a = 


g(x,y ) 0 \ „ lg(x,y ) 1 0 

0 g(x,y) r e “ \ 0 g(x,y)- 1 


r p - i (g 2 di(g 2 ) g 2 d 2 (g 2 ) 

Ag- 2 d2{g 2 ) 


r P - l/g d 2(g ) 


' k2 


g “di(g 1 ) 
2 di (g 2 ) 


~g 


2 \g- z d 1 (g z ) g- z d 2 (g z ) 


(7) 


The spin connections are calculated as 
W 21 = ~ d 2(g), 

f 

W 12 = ~ d 2 (g), 
g 

(j)“ h - 0. (other components) 


( 8 ) 


By substituting these into Eq. (|2|, we arrive at the massless 
Dirac equation on a curved surface in terms of isothermal co¬ 
ordinates: 



FIG. 1: (Color online) A helicoid in three-dimensional Euclidean 
space. Shown on the surface are the isothermal coordinates, (u, v). 


r° ^ 

W d 0 


_ F I 

y/8^-1 \ Ag^-j ’ 


(9) 


where d = d x - id v , I) = d x + id y and <// + , i// are two compo¬ 
nents of the spinor 13 . It is easy to see that Eq. |9]) has a chiral 
symmetry® since the equation has a solution with eigenenergy 
-E accompanied with a solution with eigenenergy E: 


The massless Dirac equation (|9]i on the helicoid then reads 

i / 0 d u - id v + U + \ U + \ 

a cosh u \<9„ + id v + 0 )\P-) 

(13) 

Since v is a cyclic coordinate, representing the translational 
symmetry along the axis of the helicoid, we can set 


(0 


i 


±<9 


o jy-xfg'i'-J 



( 10 ) 


Moreover, the system has a zero-mode, where and 

y[gi/j- are holomorphic and anti-holomorphic, respectively. 


III. DIRAC EQUATION ON THE HELICOID 


The helicoid, shown in Fig. [I] can be parametrized in R 3 9 
(x,y,z) as 


'x 


'sinh u cos v' 

y 

= a 

sinh u sin v 

dC 




i//±(u,v) = exp (i£v)(/)(u) ± , (14) 


in which t is physically the “angular momentum” along the z 
axis. Since the system is helical and the coordinate v spirals 
around the axis of the helicoid, ( does not necessarily have 
to be integer-valued, or cannot be quantized. This angular 
momentum, can also be regarded as the momentum along 
the helical axis. This leaves us with a differential equation for 
<p(u) below: 


0 


cosh u \d„ - £ + 


d u + {- 


tanh u 


tanh u \ ( ± 

\ / r' H 


0 


= E 


(15) 


We took a unit in which a — 1 here. We are going to stick with 
this unit system unless otherwise stated. 
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A. Solving the Equation with fixed C 


The massless Dirac equation on the helicoid, ( fl5| , is similar 
to the flat-space Dirac equation in that it takes an off-diagonal 
form. 


0 

iD-e 
cosh w 


iD( 

coshw 

0 


where 


D( — d u + £ + 


= E 


tanh u 


We square the Dirac Hamiltonian, 
of second-order differential equations: 
1 1 


0 

iD-e 
cosh w 


iD[ 
cosh w 

0 


cosh u 
1 


-D r 


-D 


= E 2 


cosh u 

D- t D,(p_ = E 2 4>_. 


coshM coshw 


(16) 

(17) 
to get a set 

(18) 
(19) 


Although it is possible to numerically solve these equations 
directly, one alternative is to cast them into a more convenient, 
Schrodinger-type form with the transformations below: 

s = sinh(M) and = <p ± (u) Vcosh u. (20) 

This transformation is chosen for geometric reasons: (i) ds 
gives the infinitesimal canonical distance in R 3 . (ii) The inte¬ 


gral measure, cosh u du, associated with the metric, Eq. (12 1 , 
is then factored out. Following these, Eq. © is transformed 
into 


d 2 

* 


Cs 


1 + S 2 (1 + V 2 )3/ 2 


m = E-ilts). (21) 


We concentrate only on the i/y _(s) = i//(,v) component of the 
spinor here and hereafter: we only have to change l into -i 
to obtain the component of the spinor. Eq. <EU- is just 
the one-dimensional Schrodinger equation with eigenenergy 
E 2 for a particle with mass m =1/2 travelling through the 
potential barrier 


Vc(s) = 


Is 


1 + S 2 (1 + S 2 ) 3 / 2 


( 22 ) 


(see Fig. |5|. Note that E : on the right-hand side of Eq. (21 1 
ensures we can concentrate only on non-negative eigenvalues 
of this eigenvalue equation. 

The potential barrier, Vr(s), approaches zero as s —> +oo 
(Fig. [2]). This implies we have no bound states with en¬ 
ergy above zero. Therefore all the states we need to con¬ 
sider (with eigenenergy E 2 nonnegative) are scattering states, 
and we have to solve the equation with boundary conditions 
in which eigenstates behave asymptotically like plane waves 
as s -> ±oo. Examples of numerically obtained wavefunc- 
tions are shown in Fig. [3] along with their asymptotes. We 
have employed the Runge-Kutta method and solved the dif¬ 
ferential equation backwards from so = 45, with its initial 



FIG. 2: (Color online) (a) The potential barrier in the transformed 
one-dimensional problem as functions of s = sinh(M) with different 
values of L (b) 2D color map of the potential against C and s. 


value conditions chosen to match the real/imaginary parts 
of i^(so) = expO'Eso) and i/j'(sq) - iE exp(iEso). We can 
see from the figure that the wavefunctions deviate from the 
asymptotic plane waves in a region very close to the axis of 
the helicoid (s ~ 0). 

We then calculate the reflection and transmission coeffi¬ 
cients, R = Re 1 '"' and T = Te !&T , using the obtained wave- 
functions and their asymptotes (see Fig. |4j. The results are 
plotted against E in Fig. [5] 

We can see from the figure that the transmission/reflection 
probabilities resemble those of a finite potential barrier, except 
that the resonance phenomena does not show up in the present 
system. This is because the potential barrier Ve(s) is some¬ 
what rounded compared with finite hard walls. The wave¬ 
length of electrons are much shorter than the typical scale at 
which the height of the potential changes. 


B. Solving the Equation with fixed E 


We also plot R and T against l (Fig. [6|. We can see from 
Fig. [6] that |7) 2 , \R\ 2 and arg T are invariant while arg R in¬ 
creases by 77 under the inversion of the angular momentum 
(t i-> -f). The reason why there is such a symmetry will be 
discussed in the next section. 
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FIG. 3: (Color online) The wavefunction with C = 0.5 at E = 0.35 
(the red, thick curves). The green, thin curves on the right side (25 < 
s < 50) represent positive asymptotes, i.e., cos (Es) and sin(£i) in 
the upper and lower panel, respectively. The blue, thin curves on 
the left side (-50 < s < -25) are fitted curves representing negative 
asymptotes. 



FIG. 5: (Color online) The transmission/reflection probabilities with 
C = 0.5 plotted against energy. 

Physical quantities tend to change significantly where the 
eigenenergy E 1 becomes comparable with the peak height of 
the potential: 

E 2 ~ ^ [|£|(£ 2 + 3) 3/2 - e 2 (i 2 - 9)] ({> (9(1)). (23) 

This corresponds to the line where physical quantities change 
considerably in Fig. [7](the green regions in (a), the purple and 
aqua regions in (b) and the purple regions in (c)). 

One interesting question worth asking here is how the trans¬ 
mission/reflection coefficients associated with some positive 
angular momentum £ (i.e., electrons spiralling upwards in the 
+Z direction) are related to those associated with the inverted 
angular momentum, -l (electrons spiralling downwards in the 
—z direction). In the transformed one-dimensional problem, 
the relation V[(s) = V-((-s) implies the following: inverting 
the angular momentum (£ i—> -€) amounts to inverting s with l 
fixed (s i-> — ,v). We are, therefore, going to discuss the change 
in physical quantities in terms of the space inversion. 

We define the transfer matrix as described in Fig. [8] 



The relations between the variables above and those in Fig. [4] 
are 1/p = R = Re' 6 * and q/p = T = Te' 6 ’. The probability 
conservation requires \p\ 2 - \q\ 2 — 1 as well as the existence of 
A such that r = q*e ,A and s = p*e' A holcP^. 

The space-inverted transfer matrix T' satisfies, by defini¬ 
tion, the following: 


FIG. 4: The boundary condition of the system as a transformed one¬ 
dimensional problem. 


C. Scattering amplitudes and phase shifts against (. and E 

Let us summarize the above results into two-dimensional 
plots against £ and E of |T'| 2 , |7?| 2 , arg T, and arg R. These are 
shown in Fig. [7] 



This means that the transmission/reflection coefficients, 7~ and 
R , of the inverted system are e‘ A /p and ~(q*/p)e ,A , respec¬ 
tively. The transmission/reflection probabilities, |7~| 2 and 7?| 2 , 
are therefore invariant under the space inversion. We can also 
show 

1 

A s; 0, arg <7 oc - -> 


0 (E —» oo) 


(26) 
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FIG. 7: (Color online) (a) 2D colour map of the reflection probability 
(|£| 2 = 1 - |7f) on l and E. (b) 2D colour map of the reflection phase 
shifts on £ and E. (c) 2D colour map of the transmission phase shifts 
on £ and E. 


FIG. 6: (Color online) The transmission and the reflection probabil¬ 
ities and the phase shifts plotted against £ at (a) E = 0.5 and (b) 
£=1.5. 


using the first Born approximation. This is shown in Appendix 
|A| This relation holds to a good approximation, according to 
our numerical calculation. This is the reason why, as we can 
see in Fig. [6] the phase shift of the transmitted wave is sym¬ 
metric about ( — 0 while that of the reflected wave is symmet¬ 
ric otherwise a jump by n at ( — 0. 

It might first seem strange that inverting angular momentum 
should bring about changes in any physical quantities. How¬ 
ever, as we are dealing with spinors liereTd only when the 
angular momentum and the spin directions are both inverted 
will the system go back to its original state. This operation 
corresponds to inverting £ and s simultaneously. 



FIG. 8: Incoming and outgoing waves. 
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IV. LOCAL DENSITY OF STATES 


We now turn to the local density of states (LDoS) of the 
system. This is one way of visualising wavefunctions when 
dealing with continuous spectra. LDoS, p(r, E), is defined, in 
terms of the Green’s function, as 

p(r, E) = - ImTrGlr,r, E). (27) 

7T 

where G(r, r', E) is the Green’s function and is a 2 x 2 matrix: 


G(r,r',E) 


(G ++ (r,r',E) G + '(r,r',£)\ 
yG~ + (r, r',E) G-(r,r',E)J- 


The Green’s function for the Dirac equation 



D + \K 


0 jw, 



(28) 


(29) 


satisfies the following relation : 

7 0 D + (r)\_ l/G ++ (r,r',£) G + -(r,r',£)\ 

\D~( r) 0 ) J\G- + (r,r',£) G-(r,r',E)J 

= —S(r - r')l, (30) 


This gives a spectral representation of the Green’s function: 


G ++ (r,r',E) = -limY 

6 -> 

n 

G (r, r',E) = -limY 


Er n i r')^(r) 

E 2 ~ El + iS ’ 
E 2 - El + iS ’ 


(31) 

(32) 


where we have assumed the completeness of {t/G(r)} and of 
{i// n (r)}, respectively. LDoS is then given by 

p(r, £) = 2 E6(E 2 - El) [\^(r)\ 2 + W n (r)\ 2 \, (33) 


or 

P(r, E) = i 2 f|iA£ 1 (r)| 2 + |^' 2 (r)| 2 ], (34) 

a=± 

if we have a continuous spectrum. Here we let 
{(/^(r), i/f*' 2 (r)} form an orthonormal basis of the eigenspace 
of D ± (r)D T (r) with eigenvalue E 2 . 

We can Fourier transform Green’s function on the heli¬ 
coid with respect to v to decompose LDoS into partial LDoS, 
Pe(s,E): 


p{s,E) = ^ P rG,£) = [|iA“;>)| 2 + \r ( p)\ 2 ] ■ (35) 

{ ( a=± 

We calculate partial LDoS numerically. As the original 
wavefunctions (asymtoptically cos (ks) and sin (ks) as s —> oo) 
do not form an orthonormal basis, we have to prepare i/'|' 1 (s) 
and tA!' 2 (s) which are orthogonal and normalised. This or- 
thogonalization of the actual wavefunctions was done using 
the Gram-Schmidt algorithm. The numerically calculated 


C = 0.5, E = 0.3 



FIG. 9: (Color online) (a) Local density of states on the helicoid on s 
with ( = 0.5 at E = 0.3. Note the symmetry with respect to the space 
inversion, sh —s. (b) Surface color plot of LDoS on the helicoid. 


LDoS is shown in Fig. [9] We can see that LDoS oscillates 
with large amplitudes in the vicinity of the axis of the heli¬ 
coid. This is natural because the potential is peaked around 
s = sinh(w) 0 (see Fig. [2]). We also plot LDoS against E 
and .s' or against f: and s. These are shown in Fig. [10] Again 
LDoS oscillates with large amplitudes around the axis of the 
helicoid in a consistent way: faster oscillations for larger E, 
and smaller amplitudes for larger {. 


V. MASSIVE CASE 


Having studied the case of massless Dirac fermions so far, 
let us move on to the case of massive Dirac fermions on the 
helicoid. One might expect that bound states could appear in 
the massive case, for which we replace E 2 with E 2 - m 2 in Eq. 
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LDoS (arb. unit) with t = 0.9 



s 


LDoS (arb. unit) at E = 1. 3 
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FIG. 10: (Color online) Local density of states against E and s (up¬ 
per) or against (. and s (lower). 




d 2 
ds 2 


i/r(s) + 


e 

1 + s 2 


Cs 

(1 + S 2 ) 3 / 2 


i/d» = (E 2 - m 2 )t/f(s). 


(36) 

We can however prove that bound states continue to be absent 
even in the massive case. 

In order to show the absence of bound states with E 2 -m 2 < 
0, we only have to find one wavefunction whose eigenenergy 
satisfies E 2 - m 2 = 0 that is everywhere nodeless. This can 
most easily be done by going back to the original Dirac equa¬ 
tion, At E = 0, this equation becomes, for d> , 


d tanh(t<) 

— + e+ —— 

du 2 


<p- = 0, 


(37) 


which can be solved analytically to give 

C j-(u 

</>-(*) = ■ (38) 

Vcosh u 


This can be translated into a wavefunction of the Schrodinger- 
type equation, (361, with its eigenenergy satisfying E 2 - nr = 
0 : 


(39) 


where s = sinh u as befor^. This means that one of the wave- 
functions of the states with E 2 - nr - 0 is nodeless (hence a 
ground state wavefunction), which means there are no bound 
states at E 2 - nr < 0. We therefore conclude the absence of 
bound states even in the case of massive Dirac particles. 

We can also show that the state with E 2 - nr = 0 is ac¬ 
tually the ground state on any surfaces if they have at least 
one cyclic coordinate. According to Eq. & the zero-energy 
wavefunctions satisfy the following differential equation: 

(d x - dy) [ y/g(x,y)i//-(x,y)] = 0 (40) 

(see Eqs. 00 for how g and <// are defined). We think of 
y as the cyclic coordinate conjugate to the momentum k, and 
the above equation becomes 

(0*+*)[VgW0-(*)] = 0, (41) 

which again can be solved analytically to give 

-lex 

4>-(x) - -=. (42) 

VsM 

In the Schrodinger-type notation the wavefunction is 

W = 0-W y[g(x) = e~ kx > 0 (43) 

in agreement with Eq. ([39} in the case of the helicoid. Hence 
the zero-energy wavefunction have no nodes, and we again 
identify it with the ground state wavefunction. 

Turning back to the massless case, the absence of (positive 
energy) bound states can be loosely accounted for the Gauss 
curvature approaching zero as u approaches positive/negative 
infinity. This is because the effective potential stemming 
from the curvature of the surface is proportional to the Gauss 
curvatursj 2 . We therefore expect the potential barrier in the 
transformed one-dimensional problem, V((s), to be roughly 
proportional to the Gauss curvature of the surface. 


VI. CONCLUSIONS AND OUTLOOK 

We have studied the Dirac electrons on the helicoid so far. 
We saw that the scattering mainly occurs just around the he¬ 
lical axis and that the bound states are absent even when we 
make fermions massive. 

One of the candidates of the physical realization of the 
present system is graphite with a screw dislocation. As our 
analysis was done using the Dirac equation and hence ignor¬ 
ing the actual lattice, the analysis is applicable when the typi¬ 
cal length scale at which the carbon honeycomb lattice bends 
is much larger than the lattice constant. 

Although we showed the absence of bound states of mass¬ 
less/massive Dirac electrons on the helicoid, we can also build 
surfaces that, unlike helicoid, allow bound states for mas¬ 
sive/massless Dirac fermions. An example of such surfaces 
is the one whose metric is given by ds 2 = (dx 2 + dy 2 )I(y 2 +A 2 ) 
with A being constant. This has a larger curvature at larger 


i//-(s) = 4>-(u) Vcosh u = e tu > 0, 
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distances and the corresponding effective potential is larger at 
greater distances. 

Further extensions of our work may include applications 
to surfaces periodic in two or three directions rather than one. 
The Dirac equation on those surfaces can no longer be reduced 
to the ordinary Schodinger-type equation unlike the case of 
the helicoid here. 

We have considered only psudo-spins in this paper. Thus 
another possible extension of the paper will be to introduce 
real spins into the analysis: this will be of great interest with 
more complicated spin-connections. This possible research 
direction is also thought-provoking as we might be able to 
consider spin current on graphene with combined pseudo- and 
real spins. 

After the submission of the present work, we came to notice 
a paper which deals with helicoidal graphene nanoribbons ®i 
While the screw axis region of helicoidal nanoribbons, which 
is the subject of the above paper, are hollow, our main in¬ 
terest in the present paper is the very region, because of its 
singularity and possible applications to graphite with a screw 
dislocation. 
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Appendix A: Derivation of Eq. { 26 } 

We take a notation in which x = s and k = E throughout this 
appendix. The one-dimensional Lippmann-Schwinger equa¬ 
tion reads 


• r*oo 

m = e ‘ kx ~k\ d y e>k ^ yW ^yy ^ A1 ) 


We can readily see that T is invariant under space inversion 17 
and that arg T ~ l/k for sufficiently large k. Similarly, 


f 


r 


dye 


,2 iky_ 


1 + y 2 


-21 *| 


dye 2iky --—t = 4ikK 0 (2\k\) ~ 4 ike 

VTT 7 


•21*1 


(A4) 


so that arg R ~ l/k for sufficiently large k. This completes the 
derivation of Eq. 


Appendix B: Nonrelativistic Fermions on the Helicoid 


We deal with with the Schrodinger equation on the heli¬ 
coid following the unpublished work by Aoki and Morimoto^ 
cited in the introduction for the purpose of comparison with 
the Dirac case. 

We first write down the Schrodinger equation on a curved 
surface with metric gif. 


h 2 Id „ d n 2 

'2m~^dq‘^ 8 dqi ~ 8^ ( *’ 


■K 2 y 


<Kq k ) 


E<p(q k ), 

(Bl) 


where (q l ,q 2 ) is the coordinate system on the surface, denoted 
hereafter as (m, v), and x 1 and k 2 are the two principal curva¬ 
tures on the surface. 

Specifically, the Schrodinger equation on the helicoid (met¬ 
ric given by becomes 


1 / d 2 

cosh 2 u \ du 2 


1 


cosh 2 u 


(j>(u) = Ec/)(u). 


(B2) 


We used the fact that the two principal curvatures of the he¬ 
licoid is k± = +l/(acosh 2 (M)) n . We also set a = 1, fr = 1 
and m = 1/2, as in the case of the Dirac particles. Here, 
rather surprisingly, the coordinate transformation for the Dirac 
equation, Eq. ( [20} can also be applied to transform this 
Schrodinger equation into a new Schrodinger equation on the 
fiat metric: 

~ ~ = ^ (B3) 


Approximation a la Born gives the first-order perturbation: 




ikx 


ikx 


2k 


2k 


dye ikVx - y W(y)e- 


.iky 


-ikx 


f 


dye 2,k} 'V(y) + e‘ 


kx r 

%J — (X 


dyViy) . 
(A2) 


Hence the first-order perturbative expression for the transmis¬ 
sion/reflection coefficients are 

• r »oo • r »oo 

T = 1-^J dyV(y),R^^J dye 2iky V (y). (A3) 


p2 

The potential of this Schrodinger equation is V(s) = ~ 


s 2 +2 


which is different from that of the Dirac case. 


4(1+j 2 ) 2 

We first note that the system does not allow for bound 
states, following a similar argument in the main body of the 
text. The absence of bound states on the helicoid, therefore, is 
just not specific to Dirac fermions. 

We also see that the potential is symmetric with respect 
to £ i—* -t and .v i-> -s. This indicates that the reflection 
phase shift does not change by n at t — 0. This is one ma¬ 
jor difference in physical quantities between the Dirac and the 
Schrodinger case. As mentioned at the end of Section IIIC 


the symmetry of the Dirac system is the combined inversion 
of the angular momentum and the (pseudo-)spin directions, 
whereas the allowed symmetry operation in the Schrodinger 
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case is to invert the angular momentum only - no spins are 
involved in the Schrodinger system. 

LDoS of the Schrodinger system, shown in Fig. 11 behaves 
also differently from that of the Dirac system. The oscillation 
amplitude of LDoS does not diminish as we go farther away 
from the helical axis unlike the Dirac case. 


FIG. 11: (Color online) Local density of states in the Schrodinger 
case with £ = 0.5 at Vfi = 0.3. 


£ = 0.5, VI = 0.3 
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